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Abstract: We analyse the steady state regime of a one dimensional Ising model under 
a tapping dynamics recently introduced by analogy with the dynamics of mechanically 
perturbed granular media. The idea that the steady state regime may be described 
by a flat measure over metastable states of fixed energy is tested by comparing vari- 
ous steady state time averaged quantities in extensive numerical simulations with the 
corresponding ensemble averages computed analytically with this fiat measure. The 
agreement between the two averages is excellent in all the cases examined, showing 
that a static approach is capable of predicting certain measurable properties of the 
steady state regime. 
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1 Introduction 



In complex systems such as granular media the energy available due to thermal fluc- 
tuations is not sufficient to cause particle rearrangement, hence in the absence of ex- 
ternal perturbations the system is trapped in a metastable state. A granular media 
may be shaken mechanically and experiments reveal a steady state regime defined by 
an asymptotic compactivity [0. The non trivial behaviour of these systems, such as 
slow relaxation dynamics and hysteresis effects arises from the fact that such systems 
have an extensive entropy of metastable or blocked states. A vertically tapped system 
of hard spheres tends to a random close packing whereas a horizontally shaken 
system crystallises 0, the stationary states obtained in these systems are not theoret- 
ically understood for the moment. Edwards has proposed that one may construct 
a thermodynamics over metastable states in the same way as Boltzmann and Gibbs 
developed over microstates, his hypothesis is that the equilibrium measure over these 
states, in gently tapped systems, is flat over all blocked states satisfying the relevant 
macroscopic constraints. Recently this scenario has been put to the test 0] where 
the authors examine a form of tapping dynamics on two systems having an extensive 
entropy of metastable or blocked configurations, the Kob- Anderson model ||^ and the 
Random Field Ising model in three dimensions. They then compared the time aver- 
age quantities obtained in the late time aging regime with those generated by a fiat 
Edwards measure numerically calculated over blocked configurations of the system. 
While the Edward's measure worked well for the Kob-Anderson model it was shown 
to be incorrect for the Random Field Ising model in three dimensions. Hence, though 
not generically true there appear to be systems where the fiat measure works. The 
developement of a thermodynamics for such systems is important to understand the 
behaviour of granular media, powders and glasses and has far reaching fundamental 
physical and practical applications. Let us add that method of simulated annealing, 
where an analogy with classical thermodynamics is used for optimization problems, is 
one of the most robust general methods of optimization used in science and industry; 
it is possible that tapping type algorithms are more efficient in certain circumstances 
(see for example [^) and thus the general theory we are searching for should clearly 
have applications well beyond physics. 

Recently tapping dynamics has been introduced on simple spin systems in order to 
draw analogies with the physics of granular media. These spin systems were spin 
glasses and ferromagnets on random thin graphs 0, |^ and the three spin model with 
ferromagnetic coupling on random hypergraphs 0. In the p-spin spherical model a 
tapping-like dynamics has been recently introduced using a time-dependent oscillating 
magnetic field |T^. Simulations on hard sphere systems with a similar tapping mecha- 
nism were also carried out in |ll|] . The general picture emerging from these studies and 



the experiment is that the compactivity of the system is increased as the tapping 
strength (the amplitude of the external perturbation) is reduced. However in the hor- 
izontal shaking experiments this is not the case. The main reason for choosing the 
spin systems studied in ||, |] is that they have an extensive entropy of metastable 
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states 0, |T2|, |T3|. That is to say that the total number of metastable states Nms is 
given by 



Nms ~ exp{NsEdw) (1) 

where SEdw is the total Edwards entropy per spin ||^. The definition of a metastable 
state depends on the dynamics of the system, here we will define a metastable state to 
be a configuration where no single spin fiip reduces the energy of the system. In this 
paper we will examine the validity of the hypothesis of Edwards for the one dimensional 
ferromagnet under the tapping dynamics of 0, §]. If one implements a zero tempera- 
ture single spin fiip dynamics such that only single spin fiips lowering the energy are 
permitted, this system has an extensive Edwards entropy and gets stuck in metastable 
states. The dynamics is made to evolve via tapping, that is to say each spin is fiipped 
in parallel with probability p G (0, 1/2]. The system is tapped after it becomes blocked 
in a metastable state and after the tap it relaxes again under the zero temperature 
or falling dynamics (by analogy with granular material in a gravitational field). Here 
the quantity corresponding to the compactivity is the energy. If one makes an analogy 
with granular media, this dynamics corresponds to a rapid relaxation dynamics be- 
tween taps, meaning that the tapping is characterised only by its strength, p, and not 
by an additional time scale. This dynamics is therefore, in a sense, the simplest case 
of realistic tapping. The Glauber dynamics of the one dimensional Ising model can 
be solved analytically (at all temperatures), however the falling dynamics defined 
here (in between taps) does not seem amenable to analytic solution. 
The tapping dynamics for this system leads after a sufficiently large number of taps to 
a steady state energy E{p) which is constant. This value of E{p) can be determined by 
a mean field theory ||] which we suspect to be exact as it reproduces certain exact 
results obtained combinatorially and also is in perfect agreement with the numerical 
simulations. Here we present numerical simulations under tapping dynamics to measure 
quantities such as correlation functions, distribution of domain lengths and energy 
fluctuations (corresponding to a specific heat) and confront the results with exact 
calculations using the Edwards measure. 

By the Edwards hypothesis we assume that the stationary measure on the tapped 
systems is 

EcOexpj-mmC]) ... 

{<^)Edw - l^j 

where {C} is the ensemble of metastable configurations and 

Z = J dENusiE) exp{-(3{p)NE) (3) 



as suggested in and recently in |T^. Here Nms{E) is the number of metastable states 
with energy per spin E and (3{p) is a Lagrange multiplier fixing the energy of the system, 
which can also be thought of as a canonical temperature coming from the tapping. 
Defining the Edwards entropy at internal energy E, SEdw{E) = \n{NMs{E))/N we find 
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that the relation determining f3{p) (the inverse Edwards temperature) is 



^ dSEdwiE) 

Pip) = \e{p) (4) 

It was shown in 0, § that E{p) is a monotonically decreasing function of p and from 



the definition @ and the calculation of SEdwiE) carried out in |]T2[ one may calculate 
the corresponding (3{p) as a function of p, which is continuous and monotonically 
decreasing. Hence p plays the role of an effective temperature which fixes the internal 
energy of the steady state. 

This does not prove that the Edwards measure is correct for this system (the hypothesis 
of Edwards does not tell us how to calculate P{p)), one now has to see if the measure 
(ID can be used to compute quantities in the steady state regime. In this paper we 
calculate using the Edwards measure: 

• The internal energy fluctuation per spin c = {{S"^) — {S)'^)/N where S is the total 
internal energy, and consequently (S) = NE. 

• The correlation functions C(r) = {SiSi-^r) and D{r) = (SiSi+iSi+rSi+r+i) ■ 

• The distribution of domain sizes. 

The results are then compared with extensive numerical simulations. 



2 Calculations using the Edwards Measure 

The Hamiltonian we consider is that for the usual one dimensional ferromagnet with 
periodic boundary conditions 

TV 

H=-Y.S^S^+l (5) 

1=1 

The partition function for the system is then given by 

N 

Z = Tr n expiPS,S,+^)9iS,^,S, + (6) 

i=l 

where the function 6{x) =0 for x < and 6{x) = 1 for x > enforces the metastability 
(each spin is stable or marginally stable in its local field). Performing a local change 
of variables to new Ising spins = Sj^j+i we obtain 

N 

Z = Tr n exp{Pai)9{ai + a,+,) (7) 

i=l 
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Hence we find that Z — Tr where T is the transfer matrix 

a 1 
1 



T-[ ? n) (8) 



with a = exp(/3). From this we find that 

E = , " (9) 

which simply determines fi. The fluctuations in the free energy are then seen to be 
given by 

c=-E{l-E'^) (10) 

The determination of the correlation function D{r) is simple as in terms of the spins 
(jj, one has D{r-) = (aiai+r)- Using the Edwards measure one obtains 

, Tr aT'^aT^-r 
D{r) = (11) 

where a is the matrix 

-( \ \) (-) 

In the large N limit one finds 

D{r) ^ {{to\a\to)f + [^J ma\t,)f (13) 

where to = (a + ^/ + 4)/2 is the maximal eigenvalue of T and ti = — 1/to the 
remaining one, \to) and \ti) denoting the respective normalised eigenvectors. Expressing 
everything in terms of the energy per spin E we find that the connected part of D{r) 
is given by 

.2. + 1 " 



D^ir) = {l-E^)[^^^) (14) 
In the same way, one can calculate the two point correlation function finding 

C{r) = C,{r) = (15) 

The matrix aT can be trivially diagonalized: aT — PDP~^, with 

with w = (a + ^/a? -4)/2. This yields 
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with 

^_ uHl-l 



Defining by P(r) the probabihty that a given domain has length r, it is easy to see 
that this is given by 

('^^1,-1) 



We find that 



1 / a \ ''"2 
p(r) = -{-] r>2 



= r < 2 (20) 

Simphfying this yields 



l + E\ f -2E ^ ''-^ 



= r < 2 (21) 

Hence, the distribution of domain sizes is geometric for n > 2, the fact that -P(l) = is a 
consequence of metastability as a domain of length 1 would be a single spin surrounded 
by two antiparallel neighbours which is unstable. 



3 Comparison with numerical simulations 

One can compare the results of numerical simulations of tapping with the above theo- 
retical ones. For a given value of the energy, let us say E, we have tapped the system 
with a strength p such that in the steady state E = E{p). The value of E{p) can be 
calculated |0, H and we recall that E{p) is maximal for p = 1/2 where it takes the value 
— l+e"^, and E{p) is monotonically decreasing for p G [0, 1/2] with limp^o+ = — 1. 
The system is tapped for a sufficiently large number of times , say tg, to ensure that the 
average of the internal energy E{t) measured becomes stationary. Once in this steady 
state regime, the quantities of interest are measured over a measurement time (number 
of taps) = 10^ . The systems were of size of 2 x 10^ spins and the results were also 
averaged over A^^ realisations. Hence, mathematically, the average value of a quantity 
A is calculated, as one would in a Monte Carlo simulation of a thermal system, as 

(^) = TrE- E Mt) (22) 

-^-^ Im t=ts+l 

In our simulations we found that for the number of sites and used here, that the 
results obtained from averaging over several systems were identical to those obtained 
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from measurements over a single system. Hence the results are in the thermodynamic 
limit and the dependence on the system size vanishes. Consequently the results pre- 
sented here are from an average over a single system of size 2 x 10^ spins. In Fig.(|l]), 
we compare the fluctuation of the energy c calculated using Edwards' measure, as a 
function of E, against those obtained from the simulations, the agreement is very good. 
For small values of p the error bars in our measurements are very small and the agree- 
ment with Eq. (|10D is excellent. For larger values of p the error bars are large as the 
statistical fluctuations are larger, however from fig.(|l]) we see that the value given by 
Eq. ([T0| ) is within the error bars. The distribution P(r) of domain sizes, is shown in 
fig.(0) and has a perfect exponential decay for r > 2. This guarantees that the result 
Eq. (^) is in perfect agreement with the simulations (as the energy E is fixed) as 
shown in Fig. (Q). 

The correlation functions such C{r) and D{r) have also been computed numerically. 
In Fig.(^) we have plotted the results in comparison with those expected from Eqs. 
(p^) and (0). Here again, the comparison is excellent (remark that the agreement 
is better for low energies, as again the statistical fluctuations due to the tapping are 
much smaller for low p than for high p). 



4 Conclusion 

We have simulated numerically tapping dynamics on a one dimensional system and 
measured fluctuations of the energy, correlation functions and distributions of domain 
sizes. The values of these quantities expected from a thermodynamics built by using 
a flat measure over blocked configurations agrees very well with our simulation data. 
The use of the flat measure therefore allows one to accurately predict the two point 
correlation function, a particular four point function D{r) and also a hierarchy of n 
point functions (see Eq. (p!9|)) corresponding to the distribution of domain sizes P(r). 
In principal further quantities could be investigated but the numerical study would be 
hampered by statistical fluctuations. The important point here is that the principal 
quantities open to experimental determination are well predicted from Edwards' mea- 
sure. A proof of the absolute validity of the use of the flat measure seems difficult, 
there is no obvious form of detailed balance in the tapping dynamics and one would 
need to show that at fixed internal energy per spin the system explores uniformly 
the metastable states of this energy. Physically this seems quite likely in the limit of 
small tapping. Here a metastable state can be viewed as a configuration of islolated 
domain walls. The first order 0{p) effect of tapping creates domain wall diffusion 
and annihilation upon the encounter of two domain walls as in the zero temperature 
Glauber dynamics of the Ising model |T3, |16|. In addition within a domain the flip- 



ping of two consecutive spins creates a domain of length 2 which then con tributes the 
the diffusion/annihilation process mentioned previously, this process is however O(p^). 
Physically, the steady state is then reached upon the equilibrium between the creation 
of small domains of this type within larger domains and the annihilation of domains 
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driven by diffusion. Hence, as the steady state regime is characterised by an average 
number of domain walls with the annihilation and creation processes mentioned above 
in equilibrium, it seems plausible the diffusion generated by the 0{p) tapping enables 
the system to explore the configurations available in a fiat manner. 
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Figure 1: The internal energy fluctuation c{E) versus E. The solid line corresponds 
to the value obtained from Eq. ([10|) and the symbols are the results obtained from 
tapping simulations made on 5000 systems of 200000 spins 
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Figure 2: Distribution of domain lengths from tapping simulations for p = .1. The 
vertical scale is logarithmic. The slope is 6 = 0.165 ± 0.001, in excellent agreement 
with that obtained from Eq. (^11) . 
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Figure 3: Slope of ln(P(r)), b{E) (as defined by Eq. (|2l|) ) with respect to the energy. 
The sohd hne correspond to from Eq. (pT]) and the symbols correspond to the same 
numerics as Fig.(|l|). 




Figure 4: Comparison between the expected C{r) and D{r) from the theoretical calcu- 
lation with the results from numerical simulations for E = —0.78 (left) and E = —0.63 
(right). The symbols are the results of the tapping experiments and the solid lines 
correspond to the one predicted by Eq.(p!7D and Eq.(|l^. 
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